本算例演示利用Fluent研究罐式燃烧器内甲烷与空气的混合燃烧过程。 算例包含的内容: 在Fluent Meshing中对燃烧室进行网格划分使用涡耗散模型模拟燃烧过程使用稳态小火焰模型模拟燃烧过程1 问题描述计算所采用的几何模型及边界示意如下图所示。 ![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247574896_cyqpvr.jpg?imageView2/0)
压缩的一次空气通过燃烧器底部的主进气口以10 m/s的速度进入燃烧室。六个涡流叶片引导进入燃烧器内的空气,并促进其与甲烷的混合以实现适当的燃烧。甲烷通过六个燃料入口以40米/秒的速度喷射进入燃烧器。当反应混合物通过燃烧器时,二次空气通过主燃烧区下游的六个二次空气入口以6 m/s的速度进入燃烧室。这有助于提高燃烧效率,并在燃烧室内壁暴露于热反应流时冷却罐壁。燃料和氧化剂以300 K的温度进入燃烧室。 2 网格生成利用Fluent Meshing中的Watertight Geometry工作流创建计算网格。 启动Fluent Meshing,启用Watertight Geometry工作流程![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247564585_hs153g.jpg?imageView2/0) 2.1 导入几何进入模型树节点Import Geometry点击按钮Advanced Options ,设置Separate Zone By为region设置Tolerance为0.1指定File Name为can_combustor.pmdb![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247576182_4lt0q3.jpg?imageView2/0) 2.2 增加局部尺寸控制进入模型树节点Add local sizing,选择Would you like to add local sizing为yes设置Size Control Type为Face Size指定Target Mesh Size为1选择边界列表项fuelinlet, inletair1及inletair2点击按钮Add Local Sizing添加尺寸控制![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247564137_snngj2.jpg?imageView2/0) 对燃料入口边界添加尺寸控制设置 Size Control Type 为 Proximity指定 Local Min Size为0.5, 指定Max Size 为2设置Cells Per Gap 为16选择列表项对象为fuelinlet点击按钮Add Local Sizing![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247574148_px5rle.jpg?imageView2/0) 对入口叶片进行尺寸控制设置 Size Control Type 为 Proximity指定 Local Min Size为0.5, 指定Max Size 为2设置Select By 为zone选择列表项origin-solid:18, origin-solid:20, origin-solid:21, origin-solid:24 及 origin-solid:25点击按钮Add Local Sizing![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247567864_yqe14k.jpg?imageView2/0) 对入口叶片进行面尺寸控制设置 Size Control Type 为 Face Size指定 Target Mesh Size为1设置Select By 为zone选择列表项origin-solid:18, origin-solid:20, origin-solid:21, origin-solid:24 及 origin-solid:25点击按钮Add Local Sizing![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247565027_1wrhjm.jpg?imageView2/0) 2.2 生成面网格进入模型树节点Generate the surface mesh设置Minimum Size 为1设置 Maximum Size 为15设置Cells Per Gap 为4点击按钮 Generate the Surface Mesh 生成面网格![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247571464_clkswy.jpg?imageView2/0)
面网格如下图所示。 ![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247581320_d6kgws.jpg?imageView2/0) 2.3 描述几何进入模型树节点Describe the geometry,如下图所示设置点击按钮Describe Geometry![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247578882_eldk9z.jpg?imageView2/0) 2.4 更新边界进入模型树节点Update the boundaries,修改边界wallvanes的边界类型为wall点击按钮 Update Boundaries![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247580508_gv0x7v.jpg?imageView2/0) 2.5 更新区域进入节点Update the regions,保持默认设置点击按钮Update Regions![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247570659_pd487g.jpg?imageView2/0) 2.6 添加边界层参数进入模型树节点Add boundary layers,采用默认设置点击按钮Add boundary layers![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247573774_y0mu4u.jpg?imageView2/0) 2.8 生成体网格进入模型树节点Generate the volume mesh指定Max Cell Length 为7.5点击按钮Generate the Volume Mesh 生成计算网格![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247567509_yjlbe7.jpg?imageView2/0)
生成的体网格如下图所示。 ![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247575218_e7awlk.jpg?imageView2/0) 点击按钮Switch to Solution 切换至求解模式3 ED模型计算燃烧3.1 General设置保持默认设置![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247566292_hggb5q.jpg?imageView2/0) 3.2 Models设置激活能量方程![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247572599_4ts7ip.jpg?imageView2/0) 采用SST k-omega湍流模型![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247578127_shnsea.jpg?imageView2/0) 设置组分输运模型![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247568246_5y1vw2.jpg?imageView2/0) 3.3 边界条件设置1、fuelinlet边界 指定fuelinlet边界入口速度40 m/s![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247566618_uzyjhl.jpg?imageView2/0) 进入Species标签页,指定组分ch4的质量分数为1 ,表示从该边界介入计算域的组分全部为甲烷![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247577759_e69udp.jpg?imageView2/0)
2、inletair1边界 指定边界inletair1的入口速度为10 m/s![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247572928_c8idet.jpg?imageView2/0) 指定该边界入口组分o2质量分数为0.23 ,表示进入的是空气![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247567018_jn1euc.jpg?imageView2/0)
3、inletair2边界 指定边界inletair1的入口速度为6 m/s![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247574534_2jof9w.jpg?imageView2/0) 指定该边界入口组分o2质量分数为0.23 ,表示进入的是空气![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247568654_vw98bg.jpg?imageView2/0)
4、outlet设置 激活选项Average Pressure Specification![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247573343_32wcom.jpg?imageView2/0) 3.4 Methods设置激活选项High Order Term Relaxation ,其他参数保持默认设置![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247565869_uo5v0j.jpg?imageView2/0) 3.5 监测物理量监测出口位置的CO3质量分数![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247580109_z3zj7t.jpg?imageView2/0) 3.6 初始化设置采用Hybrid Initialization进行初始化![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247579727_nf8lvk.jpg?imageView2/0) 3.7 计算参数设置指定Time Scale Factor为5指定Number of Iterations为500![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247569027_aaol8t.jpg?imageView2/0)
注:参数Time Scale Factor可以用于调节伪瞬态计算的时间步长,该参数值越大收敛越快,但相应的计算稳定性也越差。 ”监测得到的出口CO2质量分数变化曲线如下图所示。 ![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247565457_4x3nxz.jpg?imageView2/0) 查看进出口质量流量![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247578509_nafht6.jpg?imageView2/0) 查看系统内能量平衡![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247576528_i5uryv.jpg?imageView2/0) 3.8 计算结果y=0面上温度分布![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247569757_cq84by.jpg?imageView2/0) y=0面上co2质量分数分布![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247563728_qvnabq.jpg?imageView2/0) y=0面上o2质量分数分布![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247570993_yk0dwv.jpg?imageView2/0) 4 稳态扩散小火焰模拟前面使用涡耗散模型模拟燃烧过程。下面将使用稳态扩散小火焰模型来模拟湍流非预混燃烧过程。稳态扩散小火焰模型可以模拟由于湍流应变引起的局部化学非平衡。 在稳态扩散小火焰模型中,化学反应发生在一个局部一维的薄层区域,称为“小火焰”,湍流火焰由这些小火焰的集 合来表示。利用详细的化学反应动力学来描述燃烧过程,假设化学反应对湍流应变反应迅速,当应变弛豫到零时,化学反应趋于平衡。尽管有趋向于平衡的趋势,小火焰模型计算结果通常比涡流耗散模型或一步或两步有限速率模型能产生更精确的结果。这是因为所有的化学细节都包括在内,使得捕捉一些更快的中间反应成为可能。为了模拟湍流混合,在运行时使用概率密度函数(PDF)表作为查找表。 4.1 Models设置修改Species Model对话框中的参数选择模型Non-Premixed Combustion激活选项Steady Diffusion Flamelet选择选项Create Flamelet点击按钮**Import CHEMKIN Mechanis m…**打开机理导入对话框![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247580904_0djegv.jpg?imageView2/0) 如下图所示导入甲烷燃烧机理文件grimech30_50spec_mech.inp![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247581880_1gweux.jpg?imageView2/0) 进入Boundary标签页,指定Fuel组分ch4的质量分数为1 ,其他参数保持默认![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247571811_ecf8c3.jpg?imageView2/0) 进入Flamelet标签页,采用默认参数设置,点击按钮Calculate Flamelets 创建小火焰![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247576939_v6z0hk.jpg?imageView2/0) 进入Table标签页,采用默认参数设置,点击按钮Calculate PDF Table 创建PDF表![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247572231_4tcw4p.jpg?imageView2/0) 利用菜单File → Write → PDF... 保存PDF表4.2 边界条件设置修改边界fuelinlet的参数,指定Mean Mixture Fraction 为1![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247569389_x47471.jpg?imageView2/0) 其他边界保持默认设置注:这里设置Mean Mixture Fraction为1,表示从该边界进入到计算域中的介质全部为燃料。 ”4.3 初始化计算采用Hybrid Initialization进行初始化![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247579727_nf8lvk.jpg?imageView2/0) 4.4 迭代计算指定Time Scale Factor为5指定Number of Iterations为500![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247569027_aaol8t.jpg?imageView2/0) 4.5 计算结果y=0面上温度分布,可以看到计算得到的最高温度要比ED模型低![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247570206_uq6y5e.jpg?imageView2/0) y=0面上co2质量分数分布![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247579313_0o58jy.jpg?imageView2/0) y=0面上o2质量分数分布![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247577354_d43h4s.jpg?imageView2/0) 两种计算方法监测得到的出口位置CO2浓度比较![图片](https://public.fangzhenxiu.com/fixComment/commentContent/imgs/1624247575734_02uycw.jpg?imageView2/0)
相关文件: 链接:https://pan.baidu.com/s/1r4usHetiGIo4YSQik2DBmw 提取码:3au3 ”
|